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We have demonstrated displacement- and frequency-noise free laser interferometry (DFI) by par- 
tially implementing a recently proposed optical configuration using bi-directional Mach-Zehnder 
interferometers (MZIs). This partial implementation, the minimum necessary to be called DFI, has 
confirmed the essential feature of DFI: the combination of two MZI signals can be carried out in 
a way that cancels displacement noise of the mirrors while maintaining gravitational wave signals. 
The attained maximum displacement noise suppression was 45 dB. 
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Gravitational waves are ripples of space-time curva- 
ture. For an array of test masses separated from each 
other by a length scale ~ L, that is not too much bigger 
than the gravitational wavelength, the effect of gravita- 
tional waves can often be thought of as creating relative 
test-mass motions with order ~ Lh, or creating addi- 
tional time delays of order Lh/c for light pulses traveling 
between these test masses. One might then think that, in 
order to measure gravitational waves, we need to make 
sure that (i) our clocks must have timing errors lower 
than St ~ Lh/c, and that (ii) the test masses must follow 
geodesies within an accuracy of Sx ~ Lh. In fact, both are 
unnecessary. Point (i) was shown to be unnecessary more 
than a century ago, with the invention of interferometry 
— in which one light beam is split into two beams, each 
beam then measures the same quantity but with opposite 
sign, and then the two beams are subtracted coherently 
to eliminate the common fluctuations in light frequency 
while retaining the differential signal. An extended appli- 
cation of this idea in gravitational-wave detection is the 
so-called Time-Delay Interferometry (TDI) 1] . Point (ii) 
is more specific to the problem of gravitational-wave de- 
tection, and was addressed only recently. As was shown 
theoretically in Refs. when each test mass in an TV- 

test-mass array sends and receives light pulses from all 
other test masses, and if the array contains enough test 
masses [N > (d + 2), where d is the number of spatial 
dimensions of the array], then there exist combinations 
of time delays (i.e., signals) which do not sense test-mass 
motions or timing noises, but do (usually) sense gravi- 
tational waves. The existence of such noise-free signals 
relies on the fact that the gravitational wave contribu- 
tion to light pulse timing delays takes a form different 
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from that of test-mass motions — we therefore can com- 
bine the pulse time delays in such a way as to eliminate 
the common test-mass motions, while retaining differen- 
tial gravitational- wave signals. We will call configura- 
tions that cancel both timing and displacement noises 
Displacement- noise-Free Interferometry (DFI). 

Recently, practical optical designs of DFI using laser 
interferometry in 2- and 3-dimensions have been pro- 
posed In these configurations, the conventional, 
equal-arm Mach-Zehnder interferometer was used as a 
building block to eliminate laser noise. Four such Mach- 
Zehnder interferometers were combined, in such a way 
that they form two pairs of counter-propagating Mach- 
Zehnders. Within each pair, the two Mach-Zehnders 
share the same beamsplitters and folding mirrors; sub- 
traction of their outputs balances out displacement noise 
from motions of the folding mirrors. The two pairs share 
the same beamsplitters, which allows the elimination of 
beamsplitter displacement noise. It should be mentioned 
that DFI configurations only work for non-zero frequen- 
cies, because at nearly zero frequency gravitational waves 
are indistinguishable from relative mirror motions. For 
configurations studied so far (without optical cavities), 
the transfer function from h to effective phase shift is 
(uiLh/c)(£lL/c) a with u> the optical frequency, f2 the 
gravitational-wave frequency, and a = 2 for 3-D config- 
urations, and a — 3 for 2-D configurations. This means 
we need to observe at frequencies near c/L in order to 
have a transfer function comparable to conventional con- 
figurations. For instance, for the ground-based interfer- 
ometers with arm lengths limited to several kilometers, 
c/L becomes around 100 kHz, so it seems difficult for 
current DFI configurations to compete with them, which 
use optical cavities to increase the effective arm length 
targeting at gravitational waves at below kHz frequen- 
cies. Nevertheless, we speculate that it will eventually 
find realistic applications in future detectors, or at least 
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become one factor to consider during the invention of new 
detector configurations. Another thing which should also 
be mentioned is that DFI can cancel out frequency- and 
displacement-noise of the component optics completely, 
however, any practical noises other than these, such as 
an intensity noise of the laser, noise arising from the scat- 
tered light and electronic noise are potential noise sources 
in the DFI signal. Therefore, in practical application of 
DFI, these noises should be carefully suppressed to low 
enough levels. 

In this paper, we study DFI experimentally using a sin- 
gle pair of counter-propagating Mach-Zchndcrs, demon- 
strating the elimination of folding mirror displacement 
noise and the response to gravitational waves. We use 
lab-scale devices, with L ~ 1 meter, which indicates a 
frequency scale of ~ 100 MHz. Due to this high oper- 
ation frequency, we simulate motions of folding mirrors 
with the phase shift induced by an EOM (Electro-Optical 
Modulator) located at the position of the folding mirror. 
We place another EOM systematically at various other 
locations on the optical path, and use the phase shift 
it induces to map out the Green's Function of the de- 
vice's gravitational-wave response. We can do so because 
general relativity predicts that the effect of gravitational 
waves on our interferometers is equivalent to that of a 
medium with time- and location-dependent refractive in- 
dex, and because we do not intend to test general rela- 
tivity in this experiment. In order to see this, we write 
the phase shift imposed by a gravitational wave (in the 
so-called Transverse- Traceless, or TT gauge) on a light 
beam emitted from Xo with spatial direction N: 

4^{t) = - / N.NjhJ^it + l/c, x + Nl)dl . (1) 
c Jo 

Here t is the emission time, hJl T (t, x) is the gravitational- 
wave metric perturbation, see Eq. (6) of Ref . [4J . Except 
at places where the beams intersect (negligible) , we have 
a unique N vector, up to a sign difference (because we 
have counter-propagating beams overlapped with each 
other), at all points the light might travel through. This 
leads to a truly unique NiNj at each point, and the effect 
of the gravitational wave is the same as that of a medium 
with the following refractive-index distribution: 

n^(*,x)=2V i (x)JV J (x)^(t,x). (2) 

As depicted in Fig. [1] the full 2-D DFI configuration 
(essentially a squished version of the 3-D configuration in 
Ref. 0]) is composed of two pairs of counter-propagating 
MZIs, for a total of four MZIs. Mach-Zehnder interfer- 
ometers of equal arm length operated on the dark fringe 
are insensitive to laser frequency noise, and by super- 
imposing two counter-propagating Mach-Zchnders on the 
same optical path, the displacements of the folding mir- 
rors are sensed redundantly (once by each direction of the 
MZI) and can thus be unambiguously removed from the 
signal. An additional pair of counter-propagating MZIs 
can then be added, which shares beamsplitters with the 
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FIG. 1: (color online) Schematic of one possible 2-D opti- 
cal configuration of a displacement-noise-free interferometer. 
Light beams which follow common paths (thick, thin) travel 
in opposite directions, and the four MZIs share the two beam- 
splitters (BSI and BS2). The motion of each folding mirror 
(FM) is sensed by two MZIs, and the motion of the beam- 
splitters is sensed by all four MZIs. 



first pair-this second pair allows redundant sensing of 
beamsplitter displacements. As discussed in j4j, the sig- 
nals from the four MZIs can be combined in such a way 
as to cancel the displacements of all the optics while re- 
taining sensitivity to gravitational radiation. 

For this experiment we have constructed a partial-DFI 
composed of a single pair of counter-propagating MZIs 
to demonstrate the cancellation properties of the DFI. 
In Fig. [H the inner pair of MZIs, composed of BSI, BS2, 
FM1, FM2, has equal-length arms with folding mirrors 
located at the mid-point of each arm. As the folding mir- 
rors experience some displacement noise, the light field 
incident on those mirrors is phase modulated. Because 
the folding mirrors are at the mid-point of the arms, and 
there is light traveling in both directions in the MZI, 
these displacement noise based phase modulations will 
arrive at output 1 and output 2 simultaneously. Direct 
subtraction of the outputs 1 and 2 then allows one to can- 
cel the motion of the folding mirrors. Any signals which 
phase modulate the light and which do not arise at the 
midpoint will not arrive simultaneously at the outputs, 
and thus will not be fully cancelled. 

The transfer function of the displacement motion of 
folding mirror (FM1) to the signal port is 

Hm{a,b) = ±(u)/c)exp(-ih Q/c) (3) 
H D2 ( a ,b) = ±(w/c)exp(-i/i Q! n/c) (4) 

where c is the speed of light, to is the laser frequency 
and £1 is a fourier frequency of the displacement motion. 
The signals from output ports na and nb have opposite 
sign for phase variation, depending on the fringe condi- 
tion where the MZIs are controlled. When the folding 
mirror is located at the exact center of the MZI arm 
(suppose li a = l\p = I), all four signals display an identi- 
cal frequency dependence in their transfer function (with 
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different sign) as 

#D(i,2)(o,6) = ±(^/c) exp(-iZfi/c) (5) 

which enables the cancellation of signals due to displace- 
ment noise. 

In contrast, the frequency response to the gw signal in 
a simplest case is given as 

H gw (i,2) = ±(4u;/if2) exp(-i/fi/c) sin 2 (Zfi/2c) (6) 

for symmetric MZIs (l na — l n p — I), supposing normal 
incidence of gws with most effective polarization. As a 
substitute for exact simulation of the gw signals, one can 
generate a phase modulation signal at a location that is 
not the mid-point of an arm, and thus map the response 
of the interferometer to non-symmetrically located sig- 
nals, by observing the difference in arrival times between 
the two detection ports. 

The response function to such a signal is as given in 
Eq.Q. So, as a result, the signal combination of MZI 1 
and MZI 2 is given as 

H sgw = ±(w/c)(exp(-iZ la! Q/c) - exp(~il 10 n/c)) (7) 

which shows complete cancellation of the displacement 
motion of the folding mirror, while allowing simulated- 
gw signals to remain. 

The practical experimental setup is shown in Fig. [2] 
The laser source is a commercial solid-state Nd:YAG laser 
(Lightwave Model-126) yielding 100 mW at 1064 ran. 
The output beam is split into two, and each beam in- 
troduced to two MZIs after passing through Faraday iso- 
lators (FI), which allow detection of light exiting the 
counter-propagating MZI. The light paths of both MZIs 
are carefully adjusted and superposed on each other so 
that the two MZIs can exactly share the displacement 
noise and the effect of simulated-gw signals. The arm 
length of the MZIs affects the frequency where the MZI 
response to gw signal is maximized; in this experiment, 
the arm lengths were chosen to be around 3.6 m, to 
lower the peak frequency to several tens of MHz. To 
enable high frequency operation, the phase modulation 
due to both displacement noise and a gw signal were 
simulated using electro-optic phase modulators located 
on one of the MZI arms. The displacement noise simu- 
lator (EOM1) should be located at the exact center of 
the arm, so the position of EOM1 was carefully tuned, 
while the gw simulator (EOM2) was placed roughly 11cm 
away from the BSl. The location of the EOM2 also de- 
termines the frequency of maximum sensitivity of DFI to 
the simulated-gw signal: a more asymmetrically placed 
EOM will give a lower peak frequency. 

One of the output ports of MZI 1 (output la) was mon- 
itored with a DC-detector (DCPD) and an error signal 
created by subtracting a static offset; this signal was fed 
back to a PZT-actuated folding mirror (FM2) after ap- 
propriate filtering to give a mid-fringe locking control. 
The control bandwidth was very low, around 400 Hz, 
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FIG. 2: (color online) The experimental setup of partial DFI. 
MZI 1 and MZI 2 share a common laser source and Mach- 
Zehnder interferometer parts. The scale of the MZI was about 
4 m, for which purpose the optical paths of the MZI arms had 
to be folded several times on the optical table, which does not 
affect the DFI features in this experiment. 



so that the simulated signal by EOMs would not be sup- 
pressed by the fringe control feed back loop in the higher, 
more interesting frequency band. Once the fringe of the 
MZI 1 is controlled, that of MZI 2 is also automatically 
controlled because the two MZIs share common optical 
paths. 

Other output lights (output lb and 2a) were re- 
ceived with high-speed photo detectors (New Focus 1611- 
AC) to monitor differential optical path length varia- 
tions. The polarization optics (HWP: half wave plate and 
PBS: polarizing beam splitter) just in front of both high 
speed detectors act as an optical attenuator, which serves 
as a gain compensator for imbalanced outputs of the two 
detectors. The signals from high-speed photo detectors 
HPD1 and HPD2 were summed with a power combiner 
to produce DFI output signal, which was then monitored 
with a network analyzer (Anritsu MS4630B). 

The DFI features of bi-directional MZIs were demon- 
strated with a transfer function measurement from 
noise/signal simulators to DFI output. A swept-sine 
noise signal was provided by the internal oscillator of the 
network analyzer and was applied to the EOM(s). 

The results for displacement noise suppression is shown 
in Fig. [31 The transfer function of both detectors include 
the response functions of the simulating EOM, the pho- 
todetectors and optical and electronic phase delays from 
EOM through to the power combiner. Both amplitude 
and phase for the two signals were tuned to match each 
other using optical attenuators and path lengths adjust- 
ments so that the displacement noise signal disappears in 
the DFI signal. About 30 dB of suppression of displace- 
ment noise was attained in a wide frequency region, while 
the maximum attained suppression was 45 dB at a partic- 
ular frequency band, achieved when we tuned for maxi- 
mum suppression in that region. Any imbalance between 
the two MZI signals determines the suppression ratio; it 
is believed that the subtly different frequency response of 
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FIG. 3: (color online) The response functions from displace- 
ment noise to output signals. Dashed-red and dotted-blue 
curves are response to single detector HPD1 and HPD2, giv- 
ing almost identical response in amplitude while opposite in 
sign. Therefore summation of these two signals gives can- 
cellation of displacement noise. The suppression depends on 
a delicate tuning of the balance between two signals. Ap- 
proximately 30 dB of suppression can be seen in a wide fre- 
quency band from 10 kHz to 300 MHz in the green curve, 
meanwhile, the orange curve shows maximum suppression ex- 
ceeding 45 dB at around 25 MHz. 
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FIG. 4: (color online) The response functions from simulated- 
gw inputs to output signals. HPD1 and HPD2 response show 
almost the same response again in amplitude, whereas the 
phase delays in the higher frequency band are essentially dif- 
ferent from that of displacement noise, which is the very rea- 
son why the summed signal (DFI signal) remains gw-sensitive 
as is shown in the solid-green curve. 



our two photodetectors, which cannot be compensated by 
simple optical attenuators or path length tuning, was the 
limiting factor here. 

On the other hand, for the simulated-gw signal, the 
measured response function is shown in Fig. 3J Single 
detector outputs show response function in amplitude 
similar to to that for displacement noise; however, phase 
responses are different which is what enables gw signal 
detection with DFI. As the DFI signal still contains the 
response functions of the simulating EOM and photode- 
tectors, the measured transfer function was compensated 
with single detector data to remove them, then fitted 
with the theoretical function Eq.([7]) as shown in Fig. [5] 
Fitted parameters, scale factor A = 2.16 (2 in theory) 
and h a ~ l\p — 3.50 (3.55 in scale- measurement), agreed 
well with predicted values. 
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FIG. 5: (color online) The effects of the response functions 
of EOMs and HPDs on the DFI signal (solid-green curve in 
Fig. [4} were compensated by using the response function to 
single detector, shown in red dots. This corresponds to an in- 
trinsic DFI (bi-directional MZIs for this experiment) response 
function to a simulated gw signal. The solid curve is a fitted 
DFI response giving reasonable fitting parameters. 
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